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Abstract 

The method of generating functional, suggested for conventional systems by Kadanoff 
and Baym, is generalized to the case of strongly correlated systems, described by the 
Hubbard X operators. The method has been applied to the Hubbard model with 
arbitrary value U of the Coulomb on-site interaction. For the electronic Green's 
function Q constructed for Fermi-like X operators, an equation using variational 
derivatives with respect to the fluctuating fields has been derived and its multi- 
plicative form has been determined. The Green's function is characterized by two 
quantities: the self energy E and the terminal part A. For them we have derived the 
equation using variational derivatives, whose iterations generate the perturbation 
theory near the atomic limit. Corrections for the electronic self-energy E are calcu- 
lated up to the second order with respect to the parameter W/U (W width of the 
band), and a mean field type approximation was formulated, including both charge 
and spin static fluctuations. This approximation is actually equivalent to the one 
used in the method of Composite Operators, and it describes an insulator-metal 
phase transition at half filling reasonably well. 

The equations for the Bose-like Green's functions have been derived, describing 
the collective modes: the magnons and doublons. The main term in this equation 
represents variational derivatives of the electronic Green's function with respect to 
the corresponding fluctuating fields. The properties of the poles of the doublon 
Green's functions depend on electronic filling. The investigation of the special case 
n = 1 demonstrates that the doublon Green's function has a soft mode at the wave 
vector Q = (ir, 7r, . . .), indicating possible instability of the uniform paramagnetic 
phase relatively to the two sublattices charge ordering. However this instability 
should compete with an instability to antiferromagnetic ordering. 

The generating functional method with the X operators could be extended to 
the other models of strongly correlated systems. 
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1 Introduction 



The Hubbard model is one of the basic models in the theory of strongly correlated sys- 
tems. During its forty years of lifetime numerous approaches have been proposed for the 
investigation of the possible states of the system, the spectrum of its quasi-particles and 
the collective modes, the transport properties, the different types of ordered states and the 
phase transitions among them. Such long period of development of a model which could 
look simple at a first glance — since it contains only two parameters, the bare bandwidth 
W and the on-site Coulomb repulsion U — is determined by the circumstance that the 
case U > W is of main physical significance. But just in this case the theory does not con- 
tain a small parameter. Already the first researchers tried to avoid perturbative theories 
and used different non-perturbative approaches. Starting from the pioneering works of 
Hubbard jT] [2] [H], the method of decoupling of the double-time Green's functions (GF) 
was treated successfully. The works based on projecting the equations of motion for the 
basic operators come here |1] |S] |E] • The most productive application of this approach has 
been done with the method of composite operators [Z| [Hj P [HI] used widely not only for 
the Hubbard model but also for many other models ^1] of strongly correlated electronic 
systems. The method of the spectral density moments uses in essence the cut short of the 
equations of motions for the basic operators as well ^2] [HI] ■ Also the variational method 
of Gutzwiller belongs to the non-perturbative approaches [H], and made it possible to 
investigate qualitatively the behavior of a vast class of strongly correlated systems during 
the last four decades. The method of slave particles (slave bosons) represents an im- 
portant direction of investigation also [T3] [TE] [T7j [TH] . The basic operators are expressed 
through a product of conventional Fermi and Bose operators with subsequent exclusion of 
unphysical states. The suitable choice of a slave particle representation makes it possible 
to catch the physics of low energy states in the scope of the mean field approximation. 
Unfortunately there is no standard recipe for constructing such representations, and it is 
not always clear which one among the possible representations is the most adequate. 

During the last decade the method of the dynamical mean field theory (DMFT) has 
become quite popular ^H] [2H] • By means of this method it has been possible to investi- 
gate the behavior of almost all the models in the theory of strongly correlated systems 
in the region of strong and intermediate interactions. Apparently DMFT is the most ef- 
ficient method of investigation of these systems, although not exempt from some defects: 
it demands a huge amount of computations and has problems with the description of 
collective modes (see the review |21j). We do not mention here the numerical methods 
like Quantum Monte Carlo and small cluster diagonalization, because we concentrate our 
efforts on the analytical approaches. 

We want to pay attention to one of the analytical approaches where there is a pos- 
sibility to derive a consistent perturbative theory with respect to the parameter W/U. 
Definitely, such an approach corresponds to the perturbative theory near the atomic limit. 
The approach is based on the introduction of a generating functional Z[V], describing the 
interaction of the system with fluctuating fields depending on space and time. This func- 
tional corresponds to the generalization of the partition function of the system for the 
case of interactions with external fluctuating fields. For a proper choice of the V operator 
the different GFs of the system are expressed through variational derivatives with respect 
to fluctuating fields. 

At the beginning this method was developed for a weak interaction by Kadanoff and 
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Baym [22j[23] forty years ago. It could be generalized to a strongly correlated system 
when we express the Hamiltonian through some basic operators taking into account the 
correlations (for instance the Hubbard X operators) [3] instead of the conventional ones. 
The first time such an approach has been applied to the Hubbard model was in the limit 
U — > oo (with an additional small parameter 1/N, where N is the degeneracy of the 
electronic states) j21]. Afterwards this approach has been developed farther in the works 

[2nij2Slj2ZlE21- 

Recently we have provided a general framework for the generating functional approach 
(GFA) and we have applied it to a set of basic models of spin and strongly correlated 
electronic systems: Heisenberg model, Hubbard model for U — > oo, tj-model, sd-model, 
double exchange model |27j|32|. The results of these investigations have been generalized 
in the monograph |3Sj, published in Russian, and in the course of lectures delivered in an 
international school [3"^j . 

In this paper we apply the GFA to the Hubbard model with a finite Coulomb inter- 
action U. Supposing that U is large but of the order of W we express the Hamiltonian 
of the model in terms of the X operators and calculate the electronic and bosonic GFs. 
The latter describes the two types of collective modes: magnons and doublons. 

The electronic GF is a matrix with respect to the spin index a, the index a, indicating 
the Hubbard subbands, and the index v corresponding to the particle-hole representation. 
We have derived the equation in the variational derivatives with respect to fluctuating 
fields for it. Because the basic operators do not commute on c-values, the electronic 
Green's function is characterized by two functions of four-momenta: the self-energy E and 
the terminal part A. For £ and A the equations with the variational derivatives have been 
derived too, whereas it is possible to make iterations with respect to the parameter W/U. 
Just these iterative series represent the perturbative theory near the atomic limit 
We have limited ourselves to the first and second order corrections for E and extracted 
from them a mean field type Hmf part, which includes contributions depending only 
on the wave vector k, but not on the frequency. Ea/f consists of a term giving a shift 
to the Hubbard subbands and renormalizing its width. The last term was extracted 
from the second order correction H' 2 , which is an "uncutable" term (with respect to the 
hopping matrix element), while a "cutable" term E 2 describes the dynamical interaction 
with boson-type excitations. A procedure of extraction of the static part from E' 2 was 
borrowed from the Composite Operator Method (COM) [7 8j[9 10J. The main idea of 
this approach is that bosonic correlators, describing for example static fluctuations of 
charge, spin and pair, should not be calculated by some uncontrollable approximation 
(like decoupling or use of the equation of motion), but must be determined by means of 
general properties of the electronic GF [10]. 

The GFA, restricted to the mean field approximation, and the COM, restricted to 
a two-pole approximation, have a different structure for the electronic GF. In spite of 
this, the results obtained by these two methods for different properties of the Hubbard 
model turned out to be in very good agreement. In particular, such mean field GFs give 
two quasiparticle subbands with a gap between them, which vanishes for half-filling at 
some critical value U = U c , and an insulator-metal phase transition occurs. Detailed 
comparison of the mean field approximation in GFA and COM will be discussed below. 

Using the electron GFs we found, we can calculate Bose-like GFs for plasmons, 
magnons and doublons. In this paper we study only doublons - collective modes, describ- 
ing motion of double occupied states of the lattice sites. The equation for the doublon 
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GF has been derived. This equation contains variational derivatives of the electronic GF 
with respect to the corresponding fluctuating fields, coupled with charge densities. In the 
mean field approximation for the electronic S we have obtained the closed equation for 
the doublon GF. For the paramagnetic state at half filling (n = 1) the doublon GF has 
a soft mode at momentum Q — 7r = (71", 7r, . . . ). It indicates a possible instability of the 
uniform state against a charge density wave formation. When the filling deviates from 
unity {n < 1), the pole of the doublon GF has a gap U — 2/i, thus having the activative 
character. 

The content of the paper is the following. In part 2, based on the X operators 
formalism, the GFA is constructed. In part 3 it is derived the equation of motion for 
the electronic GF in the form of equation with variational derivatives. This equation is 
decoupled into two: one for the self energy and one for the terminal part. In part 4 the 
iterations of these equations with respect to the parameter W/U are implemented and the 
GF in the "Hartree-Fock approximation" is calculated. In part 5 we formulate a mean field 
approximation and compare GFA and COM approaches. A Bose-like GF for doublons is 
calculated in part 6 with the electronic GF taken in the mean field approximation. In 
part 7 we calculate the doublon susceptibility in the hydrodynamical regime. Finally in 
part 8 we discuss the obtained results and propose suggestions for further study of the 
Hubbard model. 



2 Introduction of the generating functional 

Let us consider the conventional Hubbard model for nondegenerate states. In terms of 
the Fermi operators the model Hamiltonian is 

ija i 

where Ci a (c\ a ) is the operator of annihilation (creation) of an electron on the site i with 
spin a, riia- = c\ a Ci a is the electron number on the same site with the given spin. Under 
the condition of a strong on-site Coulomb repulsion U > zt (where t is the hopping 
matrix element for the nearest neighbors and z is the coordination number) it is useful 
to express the Hamiltonian ([TJ in terms of the X operators. The operator Xf q for the 
site i describes the transitions between the four possible states p = |0), \a), \a), |2) - 
without any electron, with one electron possessing the spin projection o or —a and a pair 
of electrons, respectively. 

The X operators could be represented through the conventional Fermi operators by 
means of the relations 

A7° = 4(l-r^), 
X a ° = A r - 
X?° = n ia (l-n is ), 

X°° = (l- 7^X1-7^) 

The operators X? and Xf a describe the correlated creation of an electron and are Fermi- 
like /-operators; X? c and Xf° describe the flip of a spin on a site and the creation of a 
pair; they are Bose-like 6-operators, respectively. The remaining X's are called diagonal. 



xf 

X? 2 



(2) 
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We note that there are the hermitian-conjugate operators (Xf q y = Xf p . The sixteen 
X operators comprise thus the whole set, forming the algebra with the corresponding 
property of the product 

X?X» = 5 ap X?. (3) 

and the permutation relations of the anticommuting type for the /-operators while com- 
mutating for the 6-operators. We note that the conventional Fermi operators are expressed 
through the linear combinations of the X operators of the /-type 



X°° - aXf 2 . (4) 



These relations express the motion of the correlated electrons in the two Hubbard sub- 
bands. 

It is convenient to introduce the two-component spinors for the the /-operators: 

9{ia) = (_^ 2 ) , W\ia) = (X°°, aXf) . (5) 
Then the Hamiltonian (JTJ) is represented as TC = TCo + 7~Li, where 

n ^\Y^e a Xr + e 2 Xry (6) 
^i = EEE *l{i<r)t aia2 {ij)V a2 {3a). (7) 
Here we added to Hamiltonian (2.1) the term ^^(— fJ> — )n i(7} where /i is the chemical 

icr 

potential and h is the external magnetic field, that is why new notation appears: e a = 
—<y\ — A*> E2 = U — 2/i. In the quadratic form ((7|) & a (ia) represents the component of 
the spinor \P(icr), (a = 1,2); in addition we have introduced the matrix 

t a/ 3{ij) = 3= J J . (8) 

Note that the index a numerates the Hubbard subbands. With the help of the rule of 
multiplication Q for X operators, one can write the permutation relations of the spinor 
/-operators: 

[V(ia)®^(ja)] + = 5 ij X?°T z L (9) 
p{ia)®V{]o)] + = 5 l3 aX«\iTy)) 
where t x ,t v , t z are the Pauli matrices, and F" is a 2 x 2 matrix, composed of X operators: 

+ AT I ). (10) 



The permutation relations between /- and 6-operators have a commutator character: 

[V(i<T 1 ),X?]_ = 6 ij <T 1 &(i<T 1 )r B 



11) 
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In other cases of permutations, relations of type (jUJ) and (fTTj) give zero. 

Thus, an anticommutator of two ^-operators is expressed either through a diagonal or 
a 6-operator, but the commutator of \P- and 6-operators is naturally a ^-operator. Note 
the two relations 

{Xf )t = Xf, (12) 

which complete the algebra of the X operators. 

Let us write the equation of motion for the /-operator. For the thermodynamical time 
t (—(3 ^ r ^ (3, (3 = 1/kT) we start from the Heisenberg equation 

!F(l<7) = -[!F(l<7),W], 

which could be written in the case of the Hamiltonian ((01) (d) in the following form 

&(ia x ) = - E?\P(i(Ti) - F 1 CTl t(ii')^(i'o"i) (14) 
- Xf 1CT1 W(ii')^(i'(fi) + CTi!? t (i'<fi)t(i'i)ir 1, Xf 2 . 

Here a double-row matrix with respect to the spinor index was introduced 

E '=( £ e^lu)- < 15 > 

Here and in the following the numerical indexes indicate the four-dimensional coordi- 
nates including the site and the time r, i.e. i = (ii, n), a summation over the primed 
indexes is understood (it is a summation over the sites i and an integration over the time 
r). And finally the value 

£(ii') = 8(n - r[)t ililf Q = t(n')9, (16) 

has been introduced, representing the matrix over the spinor indexes (the last circum- 
stance has been specified by the symbol t). 

Thus the operator ^ represents the linear combination of the /-operators, with the 
bosonic 6-operators as the coefficients, and the matrixes E and t too. 

Following the method we have applied many times to different quantum models 
27j[32j[33j[?], we introduce the generating functional 

Z[V] = Tr (e- m Te~ v ) = e*, (17) 

where T is the symbol of the chronological product and the trace is taken over the whole 
set of variables of the system. 

For the Hamiltonian (jfij) - |7jl it is convenient to choose the operator V in the form 



V = v®?Xl?WyXl? + v^Xf," + <"'Xf/ CT ' (18) 
+v°?X 2 1 ? + v 2 l ?X™. 

It represents the linear combination of the whole diagonal and 6-operators with the 
single point fields v. Thus, differentiating the equation Z[V] with respect to the different 
u's, we can express the different GFs through the variational derivatives with respect 
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to the corresponding fields. For instance, for the single particle Bose-like GFs of the 
plasmons, magnons and the doublons we have the expressions: 



= - { TXr^)v = -^J^L., (19) 



2^(12) = -(TXfXr) v = - J^lp , (20) 

= -(TX*X?) V = (21) 

Here and further symbol (. . . )y = (...e~ v ), where (...) means averading over Gibbs 
ansamble with Hamiltonian 7i. Having been introduced in such a way, the GFs are func- 
tional of the fluctuating fields. Directing these fields to zero after taking the variational 
derivatives, we shall obtain the actual GFs, describing our system. The fermionic GF 
cannot be obtained by differentiation of Z[V]($) with respect to the single-point fields 
and it is necessary to determine the equation of motion for them. 

3 Equations of motion for electron Green's function 

We make use of the general equation of motion (see Appendix) and write it for the 
expression ((T^!^)), determining the electronic GF: 

^((T^iia^l^e^)) = ((T{^ 1 (i ( x 1 ),^ 2 (2a 2 )} + e^)) (22) 

+ ((T<(ia0^ 2 ( 2 a 2 )e^)) - ((T{\P ai (101) , V}-&1 2 (202) e~ F )) . 

Let us calculate now the anticommutator and the commutator of the ^-operators in (|22|). 
According to relations (jHJ) and we have: 

*i(^2)}+ = tfia {8 nn (Fi l ) aia2 + ^< a2 Xri) , (23) 

{^(i(7i), V}^ = Wp-Vfax) + v{ iai W{iax) + axvfW^xax)^. (24) 
Here W is the double-row matrix composed with the fluctuating fields: 

"7-CVV.r)- (25) 

After the substitution of expression ()14|) and the commutators in equation (|22p. the latter 
could be represented in the form: 



V\ 



- 5x2 [5 ai<72 ((TF^e- v )) + ^^((TXf^e^))] 

+ a 1 v 1 2 T x (lT^(m 1 )^(2 ( r 2 )e- v )) (26) 

+ ((TF 1 y H(n>)V(i>a 1 )V\2a 2 )e- v )) 

+ r z t(ii')((TAf lffl ^(i'a 1 )^ t ( 2 (T 2 )e- v ')) 

+ a l ((TX^(i'a l )^(2cr 2 )e~ v ))i(i'i)iT y . 
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Here the quantity 



(101,20-2) 







wps, 



,°"i CT i^-0 



'12, 



(27) 



has been introduced, which defines the zeroth-order approximation propagator of the 
electrons in the fluctuating single-point fields. This quantity is the 2x2 matrix with 
respect to spinor indexes. Expressing the mixed GFs through the variational derivatives 
of Z[V], we can represent the obtained equation as 



= - 5 12 a 1 (a 1 a 2 )Z[V] + a 1 (a 1 a[)t(ii')((m(ra' 1 )^(2a 2 )e~ v )) 

5 



(28) 



0i 



vfr x 5 iv - ir y t{w) 



5vf\ 



iT^(va 1 )^{2a 2 )e 



Here the double-row matrix is the differential operator with respect to the single point 
fluctuating fields: 

« 1 (0-1- 0-2) = I 4w,/v ^ ^ .:~ -^7- ] , (29) 



6v? 



where 



Ft 



( 8 , 5 

v L 



^00 + 5W 



+ 6 



(30) 



5v™ 1 5vf 



Also, let us note that i is the transposed matrix of t. 

As usual, we pass from the functional Z[V] to the functional $[V] using the substitu- 
tion: 



Z[V] 



(31) 



Then, the equation (J26j) results in a direct equation for the electronic GF: 

[g$(i<7, i'0-o - (a^o-^o^tV) -a 1 (0- 1 0-;)t(n')](^(i'(7 , 1 )^( 2(72 ) e -^ 

= - 5 12 [a 1 (0 1 2 )$] - 1 (5 11 ,<(T^ t (i'0 1 )r^ t (20 2 )e- F ) 



(32) 



- 01 



rt4> + ^02) (T^{i'<fi)i{i'i)iT y V\2a 2 )e- v ). 



8vf 5vl 



We see that the equation for the GF (Tfyfy^e^) contains the anomalous GF (T\P^^e~ v ) . 
Then, it is necessary to write the equation for it, too. 
Let us introduce the matrix GF: 



C(l2 



( (T#(i<7i)firt(202)e- V ) (T<P'(i0 1 )^(20 2 )e- v ') 
WT^t( 1(7l )^t( 2(T2 ) e -^ (T^(ia 1 )^(2a 2 )e- V ) 



(33) 



The underlined numerical index 1 in the left part represents the cumulative index, con- 
taining the space-time point 1 , the spin a±, the spinor index ai and one more index ui, 
accepting two values, specifying the matrix elements so that 



(34) 
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The matrix £(12) is an 8 x 8 matrix with respect to the collection of the discrete indexes. 
A matrix of such a rank appears automatically in the Hubbard model. Its arising is 
described "normal" states (without the Cooper's pairs) but also with broken symmetries 
as well. 

The set of four equations for the GFs in (J3~3j) could be written as a single matrix 
equation: 

^ov(ll') " (Am(ll') - (AY)(ii_>)] £(1/2) = (i$)(i2). (35) 
Here we introduced the operator matrix 

/ «1 (o-l^) -(?l$ai<TjT y T - m \ 

Mil) = $12 x 1 , (36) 

where each element represents the 2x2 matrix with respect to the spinor indexes, hidden 
in the Pauli matrices and the matrix a\, having the variational derivatives with respect 
to the fluctuating fields as its elements. Besides, the equation (J33j) contains the matrix 



Y ™ = {\ - t >)J- (37) 

The value L^ v represents the double-row matrix 

t-U,«\-( G^ v {ia 1 ,2a 2 ) ai5 CT - 1(T2 5i 2 r^° 2 \ 
0v{ - } ~ K-a^M^vf G^iuruxr*) ) ' { } 

where Gq V is given by the expression (|2*7jl . and Gq V by its transposition: 

G- V (ia u 2<r 2 ) = I (--£- + 5 aia2 + W?5 nn + v^r°5 a - ia2 J 5 12 . 

If replace back in equation ()35p term with by the mixture GFs one can see that the 
matrix equation (|33j) is equivalent that derived by Plakida [2E1 121] • The equation (|33j) is of 
the same type of the equation for the single particle GF, that we derived for the Hubbard 
model in the limit U = 00 j2Z] and for the Heisenberg model as well. In the above models 
the matrix A degenerated into a scalar, but now it is a matrix with respect to the discrete 
indexes a and v, likewise the other values in ()36|) . By virtue of the noted similarity 
of the equation (J33j) with the respective equations of the models considered before we 
could expect the same structure in the solutions of these equations, in particular the 
multiplicative character of the electronic GFs. Let us represent them as a product of the 
propagator L and the terminal n parts, respectively, namely: 

£(12) = L(ii')n(i'2). (39) 

The propagator part satisfies the Dyson equation 

^ 1 (l2)=L ^(l2)-S(l2). (40) 

Let us represent the equation for the self-energy part like the sum of the two terms: 

£(12) = x'(i2) + (nr)(i2), (4i) 
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which took place for the models considered before. Then, inserting (J4U)) and (|41|) in (|39|) 
and comparing with the initial equation ([35)1 . we can obtain the two equations for II and 
E': 

n(l2) = (i$)(l2) + (YL)(4'3')i(l4')II(3'2), (42) 
S'(l2) = -(FL)(4'3')i(l4') (Xoy(3'2) - S'(3'2)) . (43) 

In obtaining these equation we have taken into account the identity 

L(3'a), (44) 

which is the generalization of the well known identity expressing the differentiation of a 
GF through the differentiation of its inverse: 

— = G—G. 

Sv 5v 

The equations ()42j) and (|43|) for the terminal and self-energy parts of the GF have 
a structure analogous to the respective equations of the other models. These are the 
equations for the variational derivatives for II and E'. The contribution E' in the self- 
energy part E is not cutable through the "line of the interaction" , representing the value 
Y . The cutable part E has been already extracted in the equation ([4*T]) like the second 
contribution. 

From the set of the equations (|3"H|) - (|4"T|) it follows an important consequence, which 
could be represented in the form of the following equation for the GF C: 

C = C' + C'YC. (45) 
Here £' is determined by the two relations: 

£ = L'U, V - 1 = - E' . 
The solution of the equation ()45j) could be written as: 

£(12) = [C'- 1 -Y]- 1 {it), (46) 

where 

C'- l = IT 1 (L$ - S') . (47) 

As it follows from the definition, the value C is not cutable through the line Y, therefore 
the equation ()45j) for the GF is the Larkin's equation, expressing a GF through an irre- 
ducible part (with respect to a line of "interaction"). From this equation it follows the 
locator representation (|46J) for the electronic GF, also. 

So, this issue is a diagrammatic justification of the multiplicative representation 
for one-particle electron GF. Similar representations for one-particle GFs in other models 
of strongly correlated electron and spin systems was discussed in details in a review [HI] . 
The equations ([4*2|) and ([4*3"j) could be solved by iterations. At the first orders with respect 
to Y we obtain: 

n(l2) = i(l2)$ + (YL)(4'3')i(l4')i(3'2)$ + ... ; (48) 



(AL)(l2) = -L(l'2') i(ll')^ _1 (2'3') 
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(49) 



+ (YL)U>3>)(YL){e> 



L(2'5' 



In (}4*8|) the operator A, acting on $, brings the mean value of the diagonal and b- 
operators; a repeated action of the operator A will produce bosonic GFs of the different 
types. An action of the operator on L~ x will result in expressions composed of different 5- 
symbols. The problem is contained in the multiplication of the matrices in the equations 
f!48|) and (@HJ), taking into account that the matrix ^(12) contains derivatives, which 
should act on the corresponding values. To fulfil the matrix multiplication accounting for 
the operator character of the several factors, we rewrite the expressions (J48|) and (J49j) in 
another form: 

II(l2) = i(l2)$ + i(l4')(FL)(4'3')i(3'2)$ + ... , (50) 



E'l 



12 



-A(l4')(YL)(i>3')L?(3>2] 



(51) 



+i( i 40(rL)(4'30A(3'60(^)(6'l0^ 1 (l'2')^(2'5')^y 1 (5'2) + ... 



In these expressions all the factors are arranged in the order of the matrix multiplication, 
but we should not forget which factors the derivatives of the matrix A act on. 



4 Iteration equation for the self-energy and terminal 
part 

According to definition (jB3J), the electronic GF C takes into account the possibility of 
states with coupled electrons. In this paper we shall consider the normal system, described 
completely by the matrix element of the electronic GF £(12), namely 

g{ia l 2a 2 )=C n {xa 1 2a 2 ). (52) 

The normal GF Q can be looked in the standard multiplicative form 

G = GA. (53) 

with G obeying the Dyson equation 

G- 1 = G# - S, (54) 

and the self-energy part being a sum of two terms, uncutable £' and cutable At: 

E = £' + At. (55) 

In equations (p)3*j) - ()55|) all quantities are 2x2 matrices with respect to spinor indexes, 
with arguments of the type 2cr 2 )• 

Iterations in general equations (|42j) and (|43|) allow to get series for E' and A, deter- 
mined by equations (J52J) - (jo3|) . Calculations of these series are done in Appendix B, and 
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here we present the results within the limit of the first two orders. We have, for the zeroth 
order of A: 

'l - (n 5 ) 
(rf) / 

where 

{rf) = (XT + X?) = (40 (57) 



(56) 



is the average number of electrons on a site with spin a. The first order correction for A 
is the following 

AT(k) " ( K{k) K{k) \ (58) 

where 



A?(*) = - J>(* + «) 



+ (Gft + G^){k + q)V^(q) + (G& + G*,)(-fc - q)V 02 (q) 



(59) 



AJ(*) = J>(fc + «) 



(G? 2 + G£)(fc + 9)^(9) 



+ (G12 + G° 22 )(k + q)V°°{q) + {G° n + (?£,)(-* - q)V^(q) 



(60) 



The quantities N S9 (k), V a °(k) and £> 02 (fc) are the Fourier transforms of the bosonic 
GFs, determined by relations (|19|) - (|21j) with 4-momentum g. Here e(fe) is the Fourier 
transform of t^, which is actually the bare electron energy in the lattice. 
The contribution of the first order in £' is given by: 



-1 1 



(61) 



where 



if = X>(*0 [GUk)-GUk)\ 



(62) 



The second order correction is equal to: 



v-i/cr 
^2 



-ipiik) —tp%(k) 



(63) 



where 



^i( A; ) = EE £ ( fc + ^( fc i + 9) 



9 fcl 



G? 1 (W(A; + g)Gf 1 (A:i + g) 



(64) 



+ GUh)g a (k + q)G\ l {k l + q) + J2 GU-Wik + q)Gf' l (-k 1 - q) 
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and the quantity (f 2 (k) is given by a change of spinor indexes 1 <-> 2 in (J64j) . Here g a (k) 
is a linear combination of the matrix elements of the electronic GF: 



<f{k) = G° u (k) + G° 21 (k) - G° 2 (k) - G° 22 (k) 



(65) 



Finally we write down the second order contribution in the cutable part of E, that 
is, the expression for E^ ed = A^f 7 . Because in momentum representation t a is equal to 
e(fe)9f, with S being the 2x2 matrix determined in (JHJ), we find, according to (}59|) : 



YT Ted = \"(k)e(k) 



1 1 
-1 -1 



(66) 



where 



X°(k) = X^k) + X^k) 
e(k + q) g°(k + q)N* 9 {q) + g*(k + q)V^(q) 

+ g*(k + q)V 02 (q) 



(67) 



r(k) = Gl x {-k) + G° 22 (-k) - G° u (-k) - GU-k), (68) 

and 

g°(k) = -g°{-k). 

We see that the correction Ef" depends neither on momentum nor on frequency and de- 
termines only a shift of electron spectrum, but it depends on spin. The second order 
corrections ^' 2 {k) and E^" ed (/c) depend both on momentum and frequency. The contribu- 
tion E^ ed is determined by the interaction of electrons with bosonic excitations, while Y! 2 
is determined by electronic GFs only. 



5 Mean field approximation 

The simplest approximation of a mean field type is the Hubbard-I, which takes into 
account a term in E equal to A t. To it, one can add a first order term E' l5 not depending 
on frequency. The second order correction E' 2 depends on the frequency, however we shall 
try to extract from it a static part by the following ansatz. 

Let us consider that in E' 2 both expressions for ipl(k), (p 2 (k) include a factor e(k + q) 
in the summation over q. So it can be factorized in the nearest neighbor approximation, 
and a term proportional to e(k) can be taken out from the static part of E 2 for the cubic 
lattice. Thus in the static approximation E 2 can be approximated by the expression 

Z' 2 (k)=( P l P l)s(k). (69) 

Here pi and p 2 are some spin dependent constants. Their expressions can be explicitly 
written out, but we will not do it, because we shall try to calculate them from some 
general conditions for electronic GFs, which should be satisfied. Such conditions were 
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formulated in works by Mancini and coworkers (see general discussion in paper ^U] and 
refs. therein), where it is developed a method using linearized equation of motion for 
composite operators. The condition is demanding that the electronic GF Q12 is equal 
to zero when arguments coincide. Below we will use this idea for the determination of 
unknown parameters pi and p%. 

First we write down the self-energy part in an approximation which includes the 
Hubbard-I term, the first order correction £' (16 1|) and TI 2 in the form (J69|) . All these 
three contributions give Emf, corresponding to a mean field approximation. So we have: 



e(k). 



(70) 



l-(n*)+pi l-(n°)+pf 
-1 I J ' \ (n s ) - pi (n s ) -p a 2 

It is clear that the first term is responsible for a shift of the Hubbard subbands, and the 
second one for a renormalization of their widths. The propagator part of the GF in the 
mean field approximation is determined by a matrix equation: 

[G°{k))- X =[Gl{k)]- x -Y,° MF {k). 

We look for a solution of the form 

(A1U(k) {Al) aP {k) 



GUk) 



iu n -E?(k) ' iu n -E°(kY 
The poles E^(k) and their residues {A^jap^k) are written in the form: 



(71) 



a )u(fc) = ~ 



K 2 ) 22 (fe) = - 

(Al 2 ) n (k) = =f 
(Al 2 ) 21 (k) = T 



1 ± 

it 



r a (k) 

WW) 
r a (k) 



2Q°{k)_ 
^ + (l-(n ff )+pS)e(fc) 

2Q°(k) 
rf + «n ff > - pf)e(fc) 



> ■ 



2Q°(k) 



(72) 



El 2 (k) = R°(k) T Q a (k). (73) 

Here 

r CT (fc) = U - [1 - 2(n ff ) + pj + pf]e(fc), 

whilst expressions for R a (k) and Q a (k) will be written later. 

The electronic GF C? ' in the mean field approximation is found with the help of the 
general relation (jo^|) 

g*(k) = g <j (/c)Aq (fc), 

where Aq(/c) is given by the matrix ([56)1 . 

The electronic GF depends on parameters /i, r/ 7 , p^ and pf, which must be 

determined in a self consistent way from the equations 

= n h «) = W r ; i, r + 0), 
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and also from equation (|62j) . determining the parameter if . Parameters p\ and p£ wm be 
determined [10] from conditions which follow from the properties of X operators, namely: 

gUh r ; i, t + 0) = {aXf X»°) = Ol 

g^,T;z } T + 0) = (XfXpa} = 0j ' [ ' 

Thus a complete system of equations for all five parameters can be written in the form: 

(n a ) + (n*) = n, (75) 



if = X>(JO[^i(*0 

fe 



(76) 
(77) 
(78) 
(79) 



(we assume homogeneous states, so all averages do not depend on site index). 

From the comparison of the last two equations we find a relation between parameters 
Pi and p\: 

pl+pZ = -(l-2(n 9 )). (80) 

Therefore, the parameter p^ can be replaced in all the expressions above. Thus the 
equations (|72j) for the residues of GF are: 



(A ff 2 )u(fc) = j 

(^, 2 ) 22 (fc) = \ 



1 ± 
IT 



U 



2Q (J (k) 
E7 



= (Al 2 ) 21 (k) = T 



2Q CT (fc) 

^ + (( 7l ff )-p?)e(fe) 



> . 



2Q <T (fe) 

Expressions for R a (k) and Q a (k), determining poles, are now equal to: 
R°(k) = -a\ - ff + (1 + pi - (n*))e(k) + 



Q a (k) = ^ V U 2 + 4 [r + «n*> - pl)e(k)Y 



(82) 



After the replacement of parameter p\, the two equations (|75j) and (J79J) reduce to 
only one, which allows to find the unknown parameter p\. Taking the summation over 
frequencies in all equations (|76|) - (|79|) . we write our system in the form: 



(« ff >4(i-g-^(i-2K», 



(83) 
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where we use the definitions of the paper [8]: 



•J „ 



1 x ^ e n (k) 
~K^2Q a (k) 



2N 



th 



2T 



em 

2T 



-th 



2T 



2T 



(84) 
(85) 

(86) 
(87) 



We have to add to them equations ([55)) - (J5o"|) and equation (J75J) for chemical potential. 

The energy of the system can be found by averaging the Hamiltonian © - © over a 
Gibbs ensemble. It is quite easy to express it by means of electronic GFs: 



i<w> = $>(*) + f/(X 22 ) 

fcc a/3 



where 



(88) 



(89) 



After substituting here the expressions for the matrix elements C?^g, we find the expressions 
for the energy (TC) and the double occupation parameter D: 



N 



(H) = U(X 22 )+ 



(90) 



E 



(X 22 ) 



D = ^K)(l-^ + (/f ff ). 



(91) 



In Fig. 1 the parameters pi and P2 are plotted as functions of electron concentration 
at different U. Such results are typical for other fixed parameters of the system. For all 
different n and U the parameter pi is positive and P2 is negative. A negative solution for 
parameter p\ was not found. The behavior of p\ is rather similar to the COM1 solution 
for the parameter p in works jHJ [H| (COM1 is a name authors |H][01 gave for the solution 
with p > 0). In Fig. 2 the concentration dependence of chemical potential is given for 
two values of U. In the same figure a COM1 solution, that we found from equations of 
paper 0, is presented for two variants of density of states for the bare electron band: 
a two-dimensional square lattice with nearest-neighbor hopping and a model density of 
states of this type: 

< W/2 
> W/2 ' 

We see that COM and our GFA give similar results. The COM1 solution for the 2D- 
system and for the model density of states are quantitatively very close, and because of 
this we shall use hereafter for simplicity the model density of states ()92jl. 



Pie) 



fl. 


c 


(0, 


c 



(92) 



16 



Slightly worse is the comparison of results for rj (Fig.3), however there is a qualitative 
coincidence of COM1 and GFA calculations. The parameter of double occupation D gives 
again a satisfactory coincidence of the two approaches (Fig.4). It is useful to show the 
dependence of r] on n in a whole electron concentration interval at different values of U 
(Fig. 5). When decreasing U a part of the curve denoted by dash lines approaches to the 
abscess line, and when U — > one can see that rj — > 0, as it should be in the case of 
noninteracting electrons. 

Calculations show that when U decreases, the value of the jump of the chemical 
potential at n = 1 decreases too, and at the same critical value U c m 1.73 W becomes 
equal to zero. This corresponds to the closing up the two Hubbard subbands, and the 
insulator-metal phase transition occurs. The evolution of the density of states of the 
quasiparticle spectrum when U changes is shown for two different bare density of states: 
the model one (|92jl . Fig.6, and the semielliptic one 



Wl-(wV , \e\<W/2 



p{e) = {*Wy ±_ W ' l fc l^^, (93) 
, \e\ > W/2 

Fig. 7. In COM1 the critical value is U c ~ 1.68W [8J, which is close to our value U c ~ 
1.73 W, obtained for the density of states (|9*2"|). 

At half-filling it is easy to get an expression for the gap between the two Hubbard 
subbands with energies E x {k) and E 2 (k): 

AE = -(l-+p 1 )w+Ju 2 +(l--p l )w 2 . (94) 



2 

From here follows the critical value U c , when AE = 0. It is equal to 

U c =^/2p~ 1 W, (95) 

so that when U > U c the system is an insulator, and when U <U C a metal. 

Compare now the two approaches for the Hubbard model: GFA and COM. The mean 
field approximations in the framework of these approaches are close to each other both 
as what regards the GFs structure and physical properties of the model calculated with 
their help. In both cases the electronic GF has a two-poles structure. The COM approach 
includes only the parameter p, which has to be found from the equation £/ 12 = 0. In the 
GFA two parameters, (p\ and p 2 ), appear, determined through two equations: Qi 2 = 
and Q21 = 0. Due to this pair of equations one of these parameters can be eliminated, 
and as a result we have only one parameter, p\. 

The physical meanings of the parameters p and p\ are close. In the COM approach 
the parameter p describes the static fluctuation of charge, spin and pair. In GFA the 
parameter p\ includes traces of static charge and spin fluctuations as well. Corrections for 
the self-energy due to dynamical interaction of electrons with bosons in both approaches 
practically coincide and correspond to SCBA. 

The equations for the determination of parameters /1, (n cr ), 77, p\ in GFA and fi, (n 17 ), 
A, p in COM are rather similar, but have different solutions. In COM at fixed external 
parameters (n, U, W) one has two solutions: with p > and p < 0, while in GFA there 
is only one solution with p 1 > (the second parameter p 2 is always negative, but it does 
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not enter in the electronic GF explicitly; but it only guarantees the satisfaction of the 
two conditions Qu = and Q21 = 0, simultaneously). A remarkable conclusion follows 
from numerical calculations with different sets of external parameters. In spite of some 
difference in GFA and COM, the calculated quantities of the model are rather close to each 
other, if in COM only COM1 solutions with p > are taken into account. The two-poles 
GF of this approximation can be used farther for the calculation of corrections to the 
self-energy H(k,uj) from the dynamical fluctuations |3U] and for bosonic GFs (magnon, 
plasmon, doublon), describing these fluctuations pi] . 



6 Boson Green's functions 

The complete system of 16 X operators contains two Bose-like operators and X® 2 
(and their conjugates X° a and Xf°), which determine the two Bose-like GFs (2.15) and 
(2.16). 

They describe propagation of a spin flip (magnon) and a dyad (doublon), representing 
the two types of the Bose-like collective modes. These GFs could be represented as the 
variational derivatives of Z[V] with respect to the fluctuating fields: 

V a 12 = - _ _ , (96) 

V 02 (i2) = — (97) 

To write the equations of motion for the GFs T> aa and V 02 we need the equations of 
motion for the Bose-like operators: 

X? = -{s s - e a )Xf - <^,(ia)^t(ii')M^) + <^(i'a)t(i'i)^Mia), (98) 



Xl 2 = -{U - 2 f i)X° 1 2 + a'(r^) a; (ia / )^t(ii')^(i'a / ). (99) 

We see that in the right hand sides of these relations ^-operators occur; therefore in 
the corresponding equations of motions for the magnon and the doublon GFs the T- 
mixed product of /- and 6-operators will appear. They could be represented as the the 
variational derivative of the electronic GF with respect to the fluctuating field v acr in 
the first case and v 02 in the second. One of the important feature of the doublon GF is 
that it includes the "anomalous" electronic GF, composed of the two operators \P(ia) and 
l F(2a), while the equation for the magnon GF should include the normal electronic GF, 
composed of the operators \P(ia) and $^(20"). By itself these anomalous GFs are equal to 
zero when the fields are absent, however their derivatives with respect to the fields v acr 
and v 02 are not equal to zero and determine the contribution in the equation of motion, 
caused by the interactions of the electronic and bosonic degrees of freedom. Now our task 
is to determine the equations of motion for the magnon and doublon GFs and to obtain 
their approximate solution. This will let us determine the spectrum of the corresponding 
collective modes. In this paper we study only the doublon GF. 

Let us derive the equation of motion for the doublon GF (}2~Tj) ; to this purpose we write 
the equation of motion for the mean value of the operator X® 2 : 

l-((TXl 2 e- v )) = ((TX° 2 e- v )) - ((T{AT, V}e~ v % (100) 
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We substitute in it the expression (f9T?)) for X® 2 , and also the relation 
{Xf, V} = -v™X™ + vf X° 2 + vf{Xf - X 22 ). 
Then our initial equation could be rewritten in the form: 

(K$y\n>)((TX°?e- v )) = -<((T(Xf - X^O) 
+ ( r't(n0^^((^'(i'O(^^) Q '(^ / )e" V ')), 

where 

TO -1 N - - (J^ + ^ - 2A* + vf - v°A 5 12 . (101) 
Taking into account the definition of the electronic GF we write its last term in the form 

-Za'Tr [3 (tC 12 ) (ia',ia')] , 

where 

£^(l(7l,2(7 2 ) = -(TVafaJVpfa)) (102) 

is the anomalous component of the electronic GF. 

The mean values ((...)) of X operators are expressed through the variational derivative 
of the functional 3>[V], and we come to the final form of the equation for the generating 
functional: 

W)" 1 M-j^ = -v? (J^ - + (t£ 12 ) • (103) 

In the same way it is possible to write the equation for ((TXf°e~ v )) and reduce it to the 
form 

|| (/Cr' (.'0 = -v? - |* ) + <^ [3 (£ 21 <) V)] ■ (104) 

Differentiating now the equation (jlU3p with respect to v 20 , and the equation (jl04|) with 
respect to v® 2 , we come to the pair of conjugate equations for the doublon GF: 

{K™ ) - 1 (nOP 02 (i'2) = (1 - n x )5 12 - a'^Tr [9 {tC 12 ) (ia',ia')] . (105) 



V 02 (i2>) {K^T ( 2 ' 2 ) = (1 - m)S 12 - a'j-^ Tr [3 (C 21 t) (2a 1 , 2a')] . (106) 

Here we introduced the number of electrons on the site, n\ = nl+nl, where n" = (c\ a ci a ). 
We see that the exact equations for the doublon GF contain terms with variational deriva- 
tives of anomalous electronic GF with respect to the fields and v 20 . To obtain a close 
equation for doublon GF we have to calculate these terms by the same approximate way. 

Let us calculate the derivative of off-diagonal (with respect to the upper spinor indexes) 
electronic GFs C 12 and £ 21 . We use the multiplicative representation ()39|) . In the normal 
state we could use the expression for the variational derivative. 

5v 20 ~ 5v 20 L ( 44 )+ L ^) 5v 20 ( 1U7 ) 
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We take the inverse propagator GF L~ 1 in the approximation, when in the general ex- 
pression (j41j) the term E' is neglected, and II is taken in the zeroth order approximation. 
It is easy to obtain the relations: 



5 [L y(3Cr 3 ,4a 4 

5Ul 2 (30-3,40-4) _ 
5vf 
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03 ^3(74^23^34^, 



<5 2 $ 



Then, within the first order approximation with respect to t the Eq. (jl(J5|) is determined 
by the expression 



03<W A G°*{z2)T X Q°*t2A) 



(108) 



5 yA iT y - zr^(t^ ff " 3 )(3'4) P 02 ( 3 ' 2 ) 



After substituting this relation into the equation l)105|) . we represent equation for the 
doublon GF in the form 



(K^y 1 (11O - -M? 2 (ii')l P 02 (i' 2 ) = (1 - nOfe + ^ 02 (i 2 ) 



(109) 



Tr 



Q(tG a ')(l2)T X g*'(2l) 



(110) 



7W° 2 (i 2 ) = -Tr ^(tG ff ')(i 2 ) (zr^ 12 - ir^(t^')( 2 i 



The index I of the terminal and the self-energy part indicates the "left" form of the 
equations for T> 02 . In the same way starting from the equation (|106jl . it is possible to 
come to the "right" form of the equation for V 02 : 



V 02 (i2>) (K^Y'^-Mf^) =(l-n 1 )5 12 + P r 02 (i2 



(112) 



P r 02 (i 2 ) = Tr $KF (21)^(^^(12) 



M?(l2 



Tr 



3?G^'( 2 i) (iT y t l2 + ir y Q(tg a 't)(i2] 



To recover the symmetry of the doublon GF let us symmetrize the equations 
(|112|) making their sum. Then, the doublon GF is equal to 



[l-n)+V 02 (q) 



where 



M m {q) = -[MT{q)+Mf{q)} 



too n -(U-2 l 2)-M 02 (q) 

1 



(113) 

(114) 
and 

(115) 



V»\q) = -[Vy{q)+V»\q)]. 
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The self-energy and the terminal part are equal to 

M 02 (q) = — ~ ^2 £ (&) |tt[S G a (k) %t v \ — Tr[S G a (k - q) ir y ] 

ka ^ 

+\ J2 £ ( k ) \ £ ( k - G a (k) ir y 3 Q & {k - q)] 



ka 



+e(k)Tr[iT y 3 Q a {k) ^sG°{k - q)\ 



(116) 



A:<7 



In the same way we can calculate the doublon GF U 2 ®. It is possible to represent the 



result of the computation in the form 

*>^_ -a-n)+V 20 (q) 



V™{q) 



:ii8i 



(119) 



zo; n + (U-2fi) -M m {qY 
where the values V 20 (q) and Ai 20 (q) are expressed through V 02 (q) and .M 02 ^): 

.M 20 (g) = -A^° 2 (-g). 

Thus we see that the condition of symmetry is fulfilled 

V 20 (q) = V 02 (-q), (120) 

or D\\ = m the coordinate space, which follows directly from the definition (J97j) for 
the doublon GF. 

After the computation of the trace in the expressions ()116|) and ()117|) . we can represent 
them in the form: 

M 02 (q) = 5>(fc)[<f(*0 ~V(k - q)] (121) 

ka 

+~ e(fc) Uk - q)g°(k) £ Q%(k - q) + e(k)g°(k - q) £ 



a/3 



a/3 



^° 2 (?) = + G ^ k )] l^(-k + q) + QU-k + q)] (122) 

ka ^ 



+[GUk) + G" 22 (k)) [GU-k + q) + Q\ x {-k + q)) 
+ [G° u (-k + q) + GU-k + q)] [QUk) + f&(*)] 

+[GT 21 (-k + q) + GU-k + g)] [QUk) + GUk)] 
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Now we calculate the expression (|121|) in the mean field approximation for the electronic 
GF. Substituting here the formula (fTTj) and summing over frequencies, we write the result 
as a sum of two contributions of the first and second order with respect to t: 



where 



M m {q)=Mf{q)+Mf{q), 

f[E°{k)\-f[Em\ 



•Mi W - 2 2^ 2 Q°(k) 



(123) 



2 tT 

1 - / [E°(k)] - / [K(k - q)] 
iu v - E° n {k) - E*(k - q) • 1 1 

where (n, m = 1, 2) 

CUk, k-q) =[(^) n - (^) 22 ](fe)[ [«)n + (Oai] (fe)(A? 1 + A* 2 ) + 

+ [(^ + (fc) (A» + A* ) } (125) 

and here (V^), A°" are determined by formulas (jSTjl and (|56|) Remarkable is the fact 
that expressions M.\ and M.2 vanish at wave vector Q = (7r, it, ...). Because .M 02 (g) 
is nothing but the self-energy of a doublon, we see from Eq. (|118|1 . that at half-filling, 
when U — 2/x = 0, a doublon is a soft mode in the vicinity of the point (tt,tt, . . .). This 
observation pushes to study its dispersion law and attenuation. 

We postpone the study of doublons at arbitrary electron concentration and fix our- 
selves on the case n = 1. We are limited now to the hydrodynamical regime. 



7 Dynamical fluctuations in the hydrodynamical regime 

It is well known that collective modes in a disordered (symmetrical) phase in the hydrody- 
namical regime are ruled by the conservation laws [SB] • Thus the spin GF T> acr should be 
determined by the total spin conservation law, while the pseudospin GF V 02 is determined 
by the pseudospin conservation law [37] [3S] [SHI |40j . The three pseudospin components 

^+ = £^•^4 p- = J2^ Q ^c n , P* = X;i(ni-l) (126) 

i i i 

with Q = (tt 7i . . . ) obey permutation relations 

[P+,H} = (2ti-U)P+, [P-,H] = -(2h-U)P-, [P z ,H]=0, (127) 

from which it is clear that at half filling (n = 1) all pseudospin components are conserved. 
This leads to the diffusion form of the pseudospin (doublon) susceptibility, which is the 
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retarned doublon GF x 02 (q,ci;). According to Kubo-Mori theory, this susceptibility is 
expressed through the memory function Ai 02 (q,u>) by the relation 



. , Ml isu u ) 

10 7W 



where we introduce a notation for the static susceptibility, y^ 2 = X° 2 (.Q, 0). 

On the other hand the memory function is expressed through the irreducible retarded 
GF of pseudospin currents (see |4*T]): 

*«"(,, W ) = - £ /m(( f 2 |-^°f. (129) 

vy ' ; ^ ./ vr uu'(oj - uj' + i8) K ' 

V -00 

Here Xf 2 means the time derivative of operator X? 2 : 

iXf = [X° 2 ,H] = (U- 2 /U )X° 2 - a'{T x W) al {ia'){W) a/ {ia'). (130) 
Further to this, we consider the half-filling case, when U — 2fi = 0. Then 

((iX° 2 (t)\-iX?(Q)T r (131) 
= (((r^) a , (ia', t) {W) a , (ia', t) \ (& + t) pi (ja') (V + t x ) p , {ja'))) irr 

-<<(T*jr)«^,f)(&)rf(*^ 

Now we use the approximation of interacting modes, well known in the relaxation theory, 
by Mori [42J: the two-particle electron correlations in expression (J131)) are decomposed 
into pair correlators and then expressed through the imaginary parts of the retarded 
electron GFs. As a result we come to the following expression, determining the memory 
function: 

M 02 (q,u) = 2j2Hk)+e(k-q)} 2 J dJ J dw 1 [f - uj') - f (out)] (132) 

fc " "' 

tr{[Img*' (q -k,u'- oux)} [Im(^g T<7 ' {k, } 
u'(u — lu' + id) 

Here Q Ta is the transposed matrix Q° '. The quantity Q a {k,uj) is the retarded electron 
GF. It can be obtained from our Matsubara GFs by analytical continuation from discrete 
imaginary frequencies into real ones: iu n — > u> + i8. 

Expression ()132j) is similar to those obtained in the interacting modes approximation 
for other dynamical susceptibilities. For example, the spin susceptibility is: 

■xTbM = {Xf\X™)) q ,„ = M %'(qu,Y (133) 

U 

By similar decoupling of the irreducible GFs of the currents we obtain: 

M^(q,Lu)=4j2Hk)-e(k-q)] 2 f du' J du x [f [u x - u') - f {u x )\ (134) 
fc 

tr{ [Img a (k -q,u x - uj')} [Im(5Sg 9 Q)(k, uj] } 



x 



lu'(uj — lj' + id) 
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It is remarkable that the memory GF for the spin susceptibility vanishes at q = 0, 
while for the doublon susceptibility it vanishes at q = Q. This difference originates from 
the total spin conservation law (Fourier component of the spin density at q — 0), while 
the component of the pseudospin density is conserved at q = Q and only for half-filling. 
There is another important difference in expressions (jl32J) and ()134|) . Arguments of the 
electron GFs appear in a different way in these expressions. This reflects the fact that 
the spin collective mode is formed through excitations of a particle and a hole, while the 
pseudospin collective mode (doublone) is formed through excitations of two particles (or 
two holes). 

Consider now the hydrodynamical limit corresponding to small frequencies us and small 
wave point q (for expression (j!32)l ). In the hydrodynamical limit uj <C vq, where v is a 
characteristic electron velocity on the Fermi surface. Under these conditions from Eqs. 
(I132|) and (j!34|) the asymptotic expressions follow: 

ImM 02 ( qi u) = -D 02 p 2 , ReM 02 (q, u) = 0, (135) 



ImM™(q, w) = -D™q 2 , ReM™(q, u) = 0, (136) 
where the coefficients of spin and pseudospin stiffness are equal 

D 02 = 27iJ2( v ( k ) e ) 2 J duj x f{uj x )tr^Img^\-k,-uj 1 )Im[^G Ta \k,uj x )^, (137) 

D™ = 4nJ2i v ( k ) e ) 2 J du 1 f\u 1 )tr^Img a (k,u 1 )Im[^(k,uj 1 )Q}y (138) 

Here e is the unit wave vector, and f'(uj) is the derivative of the Fermi function. 

Expression (|138|) is valid at arbitrary U; in the case of U 3> W it is consistent with 
the result of |H] for the t J-model. 

Notice that if we use the electron GF in the mean field approximation (without at- 
tenuation of quasiparticles) both expressions ()137J) and ()138|) vanish. It is easy to show 
that if the attenuation of quasiparticles 7 obeys the condition 7 3> vq, both expressions 
become finite. In the general case expressions (|132|) and (|134|) for the memory function 
give in the hydrodynamical limit correct asymptotic values, therefore the susceptibilities 
have the diffusion form, which is 



1 , > Dq 2 . . 

-Im X (q, u) = Xq — 7~—T2i (139) 

u uj 1 + [Dq 2 ) 



where D = D/x^ 



8 Conclusions 

We have applied the GFA to investigate the Hubbard model in the X operator represen- 
tation. This means that we discussed the case of sufficiently strong electronic correlations 
U > W. We have derived the exact equation for the electronic GF in terms of the varia- 
tional derivatives with respect to the fluctuating fields coupled with the spin 
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and charge densities. The electronic GFs represent generally an 8x8 matrix with respect 
to the three discrete indexes a, a, v. In the matrix representation the equation has the 
same structure with the GF for the Hubbard model in the limit U — ► oo, for the tJ- and 
s<i-models and for the GFs of the transverse spin components in the Heisenberg model as 
well. 

The electronic GF Q has a multiplicative character in the sense that it is expressed 
by a product of two quantities, Q = GA, where G is the propagator satisfying the Dyson 
equation with the self energy E, and A is the terminal part. From the equation for Q, a 
pair of equations with variational derivatives for £ and A follow. Their iteration generates 
a power series in the parameter W/U. This corresponds to the perturbation theory close 
to the atomic limit. The iteration corrections of the first two orders allow to formulate a 
mean field approximation essentially equivalent to that of COM. 

Taking the electronic GF in the mean field approximation we derived an equation for 
the doublon GF. The properties of the poles of the doublon GF depends substantially 
on the electronic concentration n. For n < 1 there is a pole which has a real part 
U — 2fi > 0, corresponding to the activated mode with the quadratic dispersion law. For 
n — ► 1, U — 1\i — > 0. The investigation of the special case n — 1 reveals that a soft mode 
with Q = (jr, 7r, . . . ) may exist. However at Q = (tt, n, . . . ) the paramagnetic phase of 
the Hubbard model has an instability to antiferromagnetic ordering. It means that two 
possible instabilities - doublon and magnon ones - should compete, and a final result 
concurring a type of ordering at half filling demands a farther investigations. It will be a 
subject of next study. 

The other direction is to investigate magnetically ordered states. We should go out 
of the scope of the mean field approximation and take into account the second order 
correction £2, including the interaction of electrons with magnons. The preliminary 
analysis reveals that it contains a singular Kondo-like term ~ In | u — Ep | , which, as it 
has been pointed out in the works [121131], leads to a stable ferromagnetism. After the 
extraction of the relevant term in the second order correction we could write a more exact 
equation for the magnon GF and calculate the spin-wave spectrum. All of this will form 
the subject of a next paper. 
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Appendix A 

Equation of motion for an arbitrary Green's function 

Consider an average of an arbitrary T-product of the operators A Xl B 2 , C 3 , _D 4 ,... (it 
could be X operators, spin operators or other), taken in the Heiseberg representation: 

A 1 = B UT A i ^- H \ 1 = {«i,ti}, (A.l) 
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and so on. Then the following identity is valid: 

-^iTA x B 2 CzDt . . . e~ v j) = ([tA 1 B 2 C 3 D, . . . e~ v )) 

+{{T{A 1> B 2 }C 3 D, . . . e~ v )) + {(T{A X , C 3 }B 2 D 4 . . . e~ v )) + ... (A.2) 
-((T{A 1 ,V}B 2 C 3 D 4 ...e- y )). 

Here the average over the Gibbs statistical ensemble is denoted by 

((T...e- y )) = Tr{e- m T...e- y }. (A.3) 



The relation (jA.2|) represents the result of the differentiation of the initial average with 
respect to the time T\ ascribed to the operator A\. In the right hand side of the relation 
(lA.lj) A\ represents the time derivative 

A 1 = -[A 1 ,H}. (A.4) 

The curl brackets of the kind {Ai, B 2 } mean 

{A U B 2 } = Sin-r^A^BiMn), (A.5) 

where [..., ...]± is an anticommutator or a commutator depending on the kind of the A and 
B operators. The signs of the terms in the second line of Eq. (|A.2J) are determined by the 
signs of the transpositions of /-operators in the T-product from their original place to the 
second one in the energy term. In the last term of Eq. (jA.2j) the expression {Ai, V} is a 
commutator because the operator V is implied to be boson-like. Formally this term and 
the first one could be merged and 7i + V may be considered in a sense as the Hamiltonian 
of a system immersed in fluctuating fields. 

The identity (jA.2|) could be proved by differentiating the T-product expressed through 
the 8(t — t') functions and the expansion of the exponent e~ v in a series, with subsequent 
recollection of all the terms back to the exponent. This identity serves as a basis for the 
derivation of the GFs in the fluctuating fields, defined as 

(TA x B 2 ...e ) v = (( Te -v)) ' (A ' 6) 

where ((Te _l/ )) = Z[V] is the generating functional. 



Appendix B 

Expansion of self-energy and terminal part of the electronic Green's 
function 

The normal GF Q(i2) can be represented in a multiplicative form similar to that of 
for the whole GF, namely 

0(l2)=G(ll')0(l'2) (B.l) 

where bold indexes mean 1 = {i,a,a\}. Thus quantities Q, G, A and £' (determined by 
Dyson equation (JH1J) ) are 2x2 matrices with respect to the spinor index a. 
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We find the equations for A and S from the general matrix equations (J42|) and (|43|) 
by writing down the equation for the matrix element E' 11 : 



12 



(tL n )(4'3')A n (l4') - (tL 21 )(4'3')A 12 (l4') . (Lov) (s'2) - S' n (3' 2 ) 
(tL 12 )(4'3')i n (l4') - (tL 22 )(4'3')i 12 (l4')l ■ \{L^) 21 (3'2) - S /21 ( 3 ' 2 ) 



(B.2) 



V21, 



12 



(B-3) 



(iL 11 )(4'3')A 21 (l40 - (tL 21 )(4'3')A 22 (l4')j ■ [(L^) 11 (3'2) 
(tL 12 )(4'3')i 21 (l40 - (tL 22 )(4'3')i 22 (l4') . (L-^) 21 (3'2) 



S' n (3'2) 



3'2 



For the normal phase, matrix elements L 12 = L 21 = 0, however the derivatives of them 
with respect to the fields v 02 and v 20 should not vanish, therefore in equations (jB.2|) and 
(|B.3|) we must keep such derivatives. 

In the first order in t equations (|B.2j) and (|B.3j) lead respectively to an equation for 
the self-energy of the normal GF (E' n = £'). Matrix operators A 11 and A 12 include 
the derivatives with respect to fluctuating fields which act on (L^y) 11 and (Lq V ) 12 , and 
therefore we come to the expression: 



£i(l<7,2<7) = £'"(12) 



(B.4) 



-^12 
"#12 



(r z tG)(ia,2a) + (tr y tGr x )(ia,2a 
(tGf^n) - (tGl 2 )(u 



1 -1 
-1 1 



Here in the last line we have used a more concise definition of G — G a , and also made 
use of the expression for the transposed matrix 

G(i 2 ) = -G(ai). (B.5) 

The calculation of second order contribution in the uncutable part E' 2 demands much 
more efforts, but it involves nothing else then standard and straightforward iterations of 
equations ()B.2|) and ()B.3j) . We present the final result: 



srN = -^(^')(i2) 



£f*'( 2 i) £f*'( 2 i) 
-£f s '(2i) -B™'{2l) 



;b.6) 



where 



(B.7) 



Br (21) = (tG^) (21)^(12) - G^ 2 (2i)(G^)(i2) 

^r' (21) = (fcj,) (21)^(12) - G- 1 (2i)(G-;£)(i2) 

In a similar way we calculate the terminal part A of the electronic GF. Firstly we write 
down the expression for a matrix element II 11 = A from equation (|50p. This element is 
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coupled with the off-diagonal element II 21 . We have the following pair of coupled equations 
A(ia) (B.8) 

= (i n $)(l2) + [(tL n )(4'3')i n (l4') - (tL 21 )(4'3')i 12 (l4')l A(s' 2 ) 



+ 



(tL U )(4'3')A LL (l4>) - (tL 22 )(4'3>)A U (l4>) n 2i ( 3 '2) 



n 21 fi 2 ) 



(B.9) 



(A 2i $)(l2)+ (tL LL )(4'S')A 2L (l4')-(tL 2L )(4'3')A 22 (l4') A(3'2 



+ 



(tL 12 )(4'3')A 21 (l4') - (tL 22 )(4'3')A 22 (l4') n 21 (3'2). 

From here we find the zeroth order expressions for IT 11 and II 21 

A o (l<7i,20- 2 ) = 5i2(ai(0-i<7 2 )$), 



Il (101,202) = S 12 (Ti(iT y ] 



<5$ 



(B.10) 



(B.ll) 



Substituting these expressions in equations (|B.8jl and ()B.9j) leads to the first order cor- 
rection for the terminal part: 



5v( 



20 • 



Ai(i<t,2<t) ={T z iGT z )(i(y,2a){Tnlnl) c 

+ (T z iGr z )(ia,2a)(TXfX^} 

+ (tT y tGtT y )(ia,2a) {TXfXf) , 



(B.12) 



which includes bosonic GFs determined by definitions (|T9|) - (|2T| . 

Fourier transformations of expressions ()B.4jl . ()B.6)1 . (jB.lOJ) and ()B.12j) give the final 
results for the self-energy and the terminal part of the electronic GF. They are given by 
formulas (foTj) . . (|5fi|) and (fo^j). respectively. 
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Figure 1: Dependence of parameters pi and P2 on (a) electron concentration n and (b) 
Coulomb interaction U . 
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Figure 2: Chemical potential /x as a function of electron concentration n for two intervals 
(a) < n < 1 and (b) < n < 2. 
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Figure 3: Parameter r] as function of n and U. 
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Figure 4: Parameter D = (n a n a ) of double occupation depending on n and U . 
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Figure 6: Evolution of the quasiparticle density of states at half filling depending on U 
for a model density of states in the bare band. 
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E 

Figure 7: The same as in Fig. 6 but for a semielliptic density of states (jHSJ)- 
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